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The notion of a local fractional derivative (LFD) was introduced recently for functions of a single 
variable. LFD was shown to be useful in studying fractional differentiability properties of fractal 
and multifractal functions. It was demonstrated that the local Holder exponent/ dimension was 
directly related to the maximum order for which LFD existed. We have extended this definition to 
directional-LFD for functions of many variables and demonstrated its utility with the help of simple 
examples. 




I. INTRODUCTION 

Fractal and multifractal functions and the correspond- 
ing curves or surfaces are found in numerous places in 
nonlinear and noncquillibrium phenomenon. For exam- 
ple, isoscalar surfaces for advected scalars in certain tur- 
bulence problems typical Feynman ]3|,[I| and Brow- 
nian paths , attractors of some dynamical systems |t]] 
are, among many others, examples of occurence of contin- 
uous but highly irregular (nondiffercntiable) curves and 
surfaces. Velocity field of a turbulent fluid B at low vis- 
cosity is a well-known example of a multifractal function. 
Ordinary calculus is inadequate to characterize and han- 
dle such curves and surfaces. Some recent papers 
indicate a connection between fractional calculus ||l 3h 
and fractal structure [|||] or fractal processes 
However the precise nature of the connection between 
the dimension of the graph of a fractal curve and frac- 
tional differentiability properties was recognized only re- 
cently. A new notion of local fractional derivative (LFD) 
was introduced jl9],^0| to study fractional differentiability 
properties of irregular functions. An intresting observa- 
tion of this work was that the local Holder exponent/ 
dimension was related to the maximum order for which 
LFD existed. In this paper we briefly review the con- 
cept of LFD as applied to a function of one variable and 
generalize it for functions of several variables. 

An irregular function of one variable is best character- 
ized locally by a Holder exponent. We will use the fol- 
lowing general definition of the Holder exponent which 
has been used by various authors (2l]j2^] recently. The 
exponent h(y) of a function / at y is given by h such 
that there exists a polynomial P n (x) of order n, where n 
is the largest integer smaller than h, and 



\f(x)-P n (x-y)\ = 0(\z-y\ h ), 



(1) 



for x in the neighbourhood of y. This definition serves 
to classify the behavior of the function at y. 



Fractional calculus p3|-p5|| is a study which deals with 
generalization of differentiation and integration to frac- 
tional orders. There are number of ways (not neccesarily 
equivalent) of defining fractional derivatives and integra- 
tions. We use the Riemann-Liouville definition. We be- 
gin by recalling the Riemann-Liouville definition of frac- 



tional integral of a real function, which is given by JL3 15 



d q f(x) 1 f x f(y) 

[d{x-a)]i r(-q)J a (x-y) 



T dy for q<0, (2) 



where the lower limit a is some real number and of the 
fractional derivative 



d q f{x) 



1 



d" 



— d y ( 3 ) 



[d(x - a)]i T(n - q) dx n J a (x ~ y)?-'^ 1 

for n — 1 < q < 7i. Fractional derivative of a simple 
function f(x)=x p p > — 1 is given by p^Jl5[| 



d q x p 
dxi 



T(p+1) 

r(p - q + 1) 



x p q for p > —I. 



(4) 



Further the fractional derivative has the interesting prop- 
erty (see ref fjJJ]), viz, 



d q f{px) 
dxi 



q d q f{Px) 
d(f3x)i 



(5) 



which makes it suitable for the study of scaling. Note the 
nonlocal character of the fractional integral and deriva- 
tive in the equations and (||) respectively. Also it is 
clear from equation (0) that unlike in the case of integer 
derivatives the fractional derivative of a constant is not 
zero in general. These two features make the extraction 
of scaling information somewhat difficult. The problems 
were overcome by the introduction of LFD in Jl9| . In the 
following section we briefly review the notion of LFD for 
the real valued functions of real variable. In the section 
HI we generalize this definition to real valued functions 
of many variables and demonstrate it with the help of 
some simple examples. 
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II. LOCAL FRACTIONAL DIFFERENTIABILITY 

Unfortunately, as noted in section [|, fractional deriva- 
tives are not local in nature. On the other hand it is de- 
sirable and occasionally crucial to have local character in 
wide range of applications ranging from the structure of 
differentiable manifolds to various physical models. Sec- 
ondly the fractional derivative of a constant is not zero, 
consequently the magnitude of the fractional derivative 
changes with the addition of a constant. The appropriate 
new notion of fractional differentiability must bypass the 
hindrance due to these two properties. These difficulties 
were remedied by introducing the notion LFD in Jl9| as 
follows: 

Definition 1 If, for a function f : [0, 1] — > Ft, the limit 



E>«f(y) = lim 



d"(m - f(y)) 



v d(x — y) q 



(6) 



exists and is finite, then we say that the local fractional 
derivative (LFD) of order q (0 < q < 1), at x — y, exists. 

In the above definition the lower limit y is treated as 
a constant. The subtraction of f(y) corrects for the 
fact that the fractional derivative of a constant is not 
zero. Whereas the limit as x — > y is taken to remove 
the nonlocal content. Advantage of defining local frac- 
tional derivative in this manner lies in its local nature 
and hence allowing the study of pointwise behaviour of 
functions. This will be seen more clearly after the devel- 
opment of Taylor series below. 

Definition 2 We define critical order a, at y, as 
a(y) = Sup{q | all LFDs of order less than q exist at y}. 



These definitions were subsequently generalized |20| for 
q > 1 as follows. 

Definition 3 //, for a function f : [0, 1] — > Ft, the limit 



E>"f(y) = lim 



d"(f(x) £ 



TV 



/ ( "'fa) f T 
=0 r(n+l) V x 



y) n ) 



[d{x - y)]i 



(7) 



exists and is finite, where N is the largest integer for 
which N th derivative of f(x) at y exists and is finite, 
then we say that the local fractional derivative (LFD) of 
order q (N < q < N + 1), at x — y, exists. 

We consider this as the generalization of the local deriva- 
tive for order greater than one. Note that when q is a 
positive integer ordinary derivatives are recovered. The 
definition of the critical order remains the same since, 
for q < I, (||) and (Q) agree. This definition extends 
the applicability of LFD to C 1 — functions which are still 
irregular due to the nonexistence of some higher order 
derivative (i.e. belong to class C 7 , 7 > 1). For example 
the critical order of f(x) — a+bx+c\x\ J , 7 > 1, at origin, 



according to definitions || and || is 7. It is clear that the 
critical order of the C°° function is 00. 

In it was shown that the Weierstrass nowhere dif- 
ferentiable function, given by 



W A (t) = £A<' 



-2)k , 



(8) 



k=l 



where A>l,l<s<2 and t real, is continuously locally 
fractional differentiable for orders below 2 — s and not 
for orders between 2 — s and one. This implies that the 
critical order of this function is 2 — s at all points. Inter- 
esting consequence of this result is the relation between 
box dimension s || of the graph of the function and the 
critical order. This result has for the first time given a 
direct relation between the differentiability property and 
the dimension. In fact this observation was consolidated 
into a general result showing equivalence between criti- 
cal order and the local Holder exponent of any continu- 
ous function. The LFD was further shown to be useful in 
the study of pointwise behaviour of multifractal functions 
and for unmasking the singularities masked by stronger 
singularities. 

Whenever it is required to distinguish between limits 
from right and left sides we can write the definition for 
LFD in the following form. 



E> q ± f(y)= lim 



d q F N (x,y) 
± [d±(x- y)]9 



(9) 



The importance of the notion of LFD lies in the fact 
that it appears naturally in the fractional Taylor's series 
with a remainder term for a real function /, given by, 

where x — y = A > and R q (y, A) is a remainder given 
by 



1 



T(q+l)J 



dF(y,t;q,N) 
dt 



(A-t) q dt (II) 



and 



F(y,A;q,N) 



d"(f(x) 



i-^dfl— 



f M (y) 



n=0 T(ri+1) 



A") 



[dA]< 



(12) 



We note that the local fractional derivative as defined 
above (not just fractional derivative), along with the first 
N derivatives, provides an approximation of f(x) in the 
vicinity of y. We further remark that the terms on the 
RHS of eqn© are nontrivial and finite only in the case 
when q equals a, the critical order. Osier [ p3[ constructed 
a fractional Taylor series using usual (Riemann-Liouvillc) 
fractional derivatives which was applicable only to ana- 
lytic function. Further Osier's formulation involves terms 



with negative orders also and hence is not suitable for ap- 
proximating schemes. When A < 0, a similar expansion 
can be written for E) q _f(y) by replacing A by —A. 
When < q < 1 we get as a special case 

E> q f(v) 

f(x) = f(y) + —^(x-yy + Remainder (13) 

provided the RHS exists. If we set q equal to one in equa- 
tion ( |T3| ) one gets the equation of the tangent. All the 
curves passing through a point y and having same the 
tangent, form an equivalence class (which is modelled by 
a linear behavior). Analogously all the functions (curves) 
with the same critical order a and the same ID a will form 
an equivalence class modeled by the power law x a . This 
is how one may generalize the geometric interpretation 
of derivatives in terms of tangents. This observation is 
useful in the approximation of an irregular function by a 
piecewise smooth (scaling) function. One may recognize 
the introduction of such equivalence classes as a starting 
point for fractional differential geometry. 



III. GENERALIZATION TO HIGHER 
DIMENSIONAL FUNCTIONS 

The definition of the Local fractional derivative can be 
generalized for higher dimensional functions in the fol- 
lowing manner. 

Consider a function / : M n — > M. We define 

*(y,*) = /(y + vt)-/(y), v G M n , teM. (14) 

Then the directional-LFD of / at y of order q, < q < 1, 
in the direction v is given (provided it exists) by 



(15) 



where the RHS involves the usual fractional derivative of 
equation (||) . The directional-LFDs along the unit vector 
ej will be called i th partial-LFD. 

Now let us consider two examples. 
Example 1: Let 

oo 

Wa( x ) = W\(x, y) = ^2 ^ s ~ 2 ^ k sin\ k (x + y), (16) 

k=l 

with A > 1 and 1 < s < 2. Let v = (v x ,v y ) be a unit 
2- vector. Then 

W^(x + vt) = W x {x + v x t,y + v v t) 

OO 

= ^2 A (s_2)fe smA fc (x + v x t + y + v y t). 

k=l 

$(x,i) = W A (x + vt)-W A (x) 

OO 

= ^2 ^ {s ~ 2)k [sin\ k (x + v x t + y + v y t) 
-sin\ k (x + y)] 



If we choose y = 0, i.e., we examine fractional differen- 
tiability at a point on x axis and if we choose v y = 0, 
i.e., we are looking at LFD in the direction of increasing 
x then we get 



$(x,i) = ^\ [s - 2)k [sin\ k {x + v x t) - sin\ k (x) 
fe=l 



which is known jl9| to have critical order 2 — s. 

If we keep y = but keep v x and v y non-zero, i.e., if we 
examine fractional differentiability at a point on cc-axis 
but in any direction then we get 

OO 

$(x, t)=^2 X (s ~ 2)k [sin\ k (x + v x t + v y t) - sin\ k (x)} 
fe=l 

This again has critical order 2 — s, except when v x = —v y 
(in which case it is oo since the function (f> identically 
vanishes) . 

In fact it can be shown that the given function has 2 — s 
as critical order at any point and in any direction except 
the direction with v x — —v y . 
Example 2: Another example we consider is 

OO 

Wa(x) - W x (x,y) = J2^ 2)k Mn\ k (xy), (17) 
fc=i 

where A > 1 and 1 < s < 2. 

W x (x + vt) = W x (x + v x t , y + v y t) 

OO 

= ^ (s ~ 2)k sin(X k (x + v x t)(y + v y t)). 

k=l 



$(x,t) = W A (x + vt) - W A (x) 

OO 

= ^2 ^ s ~ 2)k [sin(\ k (xy + yv x t + xv y t + v x v y t 2 )) 
-sin(\ k {xy))} 



k=l 



k=l 



Now if we choose y — and v y — then the critical order 
is oo. 

If we choose y = but v x and v y non-zero then we get 

OO 

$(x,i) = ^2 ^ 2)k [sin(X k (xv y t + v x v y t 2 ))} 

k=l 



Therefore using results of the critical order along any 
other direction v, at a point on s-axis is seen to be 2 — s. 



IV. CONCLUSIONS 

The usefulness of the notion of LFD was pointed out in 
fL9| , p0f where the considerations were restricted to func- 
tions of one variable only. It allows us to quantify the 
loss of differentiability of fractal and multifractal func- 
tions. The larger the irregularity of the functions the 
smaller is the extent of differentiability and smaller is the 
value of the Holder exponent. Local Taylor series expan- 
sions provide a way to approximate irregular functions 
by a piecewise scaling functions. In the present paper we 
have demonstrated that it is possible to carry the same 
theme even in the multidimensional case. In particular, 
the Holder exponents in any direction are related to the 
critical order of the corresponding directional-LFD. We 
note that, whereas a one dimensional result is useful in 
studying irregular signals, the results here may have util- 
ity in image processing where one can characterize and 
classify singularities in the image data. We note that it is 
possible to write a multivariable fractional Taylor series 
which can be used for approximations and modelling of 
multivariable multiscaling functions. This development 
will be reported elsewhere. Further these considerations 
provide a way for formulating fractional differential ge- 
ometry for fractal surfaces. 
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